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Abstract 

For the first time the phenomenon of cellular structure coarsening 
are consistently analysed from the positions of kinetic, hydrodynamic and 
stochastodynamic theories of nonequilibrium statistical systems. Thereby 
micro-, meso- and macroscopic levels of approach are distinguished. At 
the microscopic level the cellular structure is describe by a probability 
distribution function in a phase space of cell coordinates and of cell sizes. 
A kinetic equation for the function is written and a development to a 
hydrodynamic equation of a mesoscopic cell medium is realised. It has the 
form of a diffusion-reaction equation with a negative "diffusion" coefficient 
and with a cell interface density playing the role of concentration. Its 
analysis reveals a new effect of macroscopic patterning in the cell medium: 
existence of space-correlated stochastic fluctuations of the cell interface 
density. 

PACS numbers: 02.50.-r, 05.20.-Dd, 47.53. +r, 82.20.-w 



1 Introduction 

There are many physical systems consisting of homogeneous domains separated 
by distinct boundaries, so called cellular structures: ordered domains in alloys, 
magnetic domains, gas bubbles in soap froth and in lipid monolayers, grains in 
polycrystals, etc. jl], ||, [| [|. We have here a simple physical situation (Fig. 1 
Q ) , when the system volume V being the sum of volumes of all cells, the volume 
specific energy and the interface specific energy E are constant. The internal 
energy of the system £ is completely defined by the cell interface area S. It is 
easy to see that for any two neighbouring cells the volume decrease of the smaller 
one (down to complete disappearance) and the simultaneous volume increase of 
the greater one diminish S. This permanent topological possibility to decrease 
£ = EH makes the system intrinsically unstable. The mean cell size increases 
during this statistical process (Fig. 1), named "coarsening" (large cells "eat up" 
small ones). The changes of £ related to the elementary statistical event of cell 
disappearance are so large in comparison with the energy of thermofluctuations, 
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Figure 1: Soap froth configurations at three different times during its evolution. 



that the coarsening cannot be described by ways of usual thermodynamics g . 
Still this challenging statistical problem can be elegantly solved by methods of 
physical kinetics, as shown by Lifshitz, Slyozov and Hillert (LSH) j|, @, § , if we 
consider the cellular structure as some dynamic system. To get a cell dynamics 
equation let us following LSH approximately consider each cell as an elastic 
sphere shell embedded in an elastic environment which represents the average 
effect of cell-cell interaction. For a cell of radius a we write an interface tension 
force as F = —E(l/a — K), where the term 1/a is related to the effect of interface 
surface tension of the cell itself and the term k - to the effect of interface surface 
tension of adjoining cells (the elastic environment). The quantity k is nothing 
that a self-consistent mean of cell curvature throughout the system. Believing 
that the cell interface velocity w affected by F is determined by viscous law 
w = MF, where M is the cell interface mobility, we have 0, || 

da/dt = -D(l/a- k) (1) 

with D — ME. In the coarsening process the cells move in a size space {a} 
(some phase space) with a velocity (1). We can write a continuity equation in 
this space as 

df/dt + d{fda/dt)/da = (2) 

for a cell size distribution function f(a,t) |], [?], |). The function f(a,t) is 
normalised so that n(t) = j Q fda is the mean number of cells per unit volume; 



n d [ a d fda= 1; 
Jo 



(3) 



a(t) = J Q afda/n is the mean cell size and p(t) = S/V = (f^d/2) f Q a d ~ 1 fda 
is the mean cell interface density. Here d = 2, 3 is the cellular structure dimen- 
sionality and Qd = 2tt(c?— l)/d. Condition (3) is nothing that the conservation 
law of the total volume of all cells in unit volume. A remarkable property of 
this statistical dynamic system, discovered by LSH, is the occurrence of a stable 
steady-state regime of evolution. This regime is observed if and only if the mean 
curvature n(t) take the unique "self-consistent" value n(t) = y/2/(Dt). All other 
macroscopic system variables obey also simple scaling laws:d(t) cx p~ x {t) oc i 1 / 2 , 
n(t) oc i d / 2 and the solution of (l)-(3) has a universal scaling form 

fo(a, t) cx (l/a(t)) d+1 cp( Pl a/d(t)), (4) 

where <p(z) is a peaked function, depending on d only, and p\ sa 1 [?], [§]. 



2 



There is a great number of publications devoted to the question of " improve- 
ment" of the individual cell dynamics equation (1) (we will say microscopic dy- 
namics equation) to get a good agreement of the cell size distribution fa(a,t) 
with experimental data, essentially in the metallurgic literature related to nor- 
mal grain growth |l(| . Surprisingly small attention has been paid to 
mesoscopic and macroscopic behaviours of the cellular structures as some con- 
tinuous elastic cellular medium with viscous motion law, although the existence 
of interesting self-organisation effects is almost evident here Jll], [13]. Indeed, if 
we have a good look at Fig. 1 we can see a strong size correlation of adjacent 
cells. Near a small cell the probability to find another small cell is grater then a 
large one, and vice versa, near a large cell the probability to find another large 
cell is grater then a small one. In other words we observe here space-correlated 
fluctuations of the mesoscopic cell interface density /?(x, t) . Qualitatively this 
effect can be easy analysed. Let us consider a nonhomogeneous cellular struc- 
ture with the macroscopic mean interface density p(t) and the mesoscopic one 
p(x., t) = p(i) + 5p(x, t) , weakly varying in space. To be more precise, let us as- 
sume that Sp(K, t) << p(t) and the wavelength A of spatial variations of p(x, t) 
is much greater than a w l/p(t). Let us cut out a medium volume Vl having the 
size L << A, but still including many cells. Each element ds of the surface E^, 
enclosing Vl, is affected by the interface tension of the cells, surrounding Vl- 
The total interface tension force acting on Vl is § Epds = J EWpdV. In a ho- 
mogeneous cell medium this force equals zero. In the nonhomogeneous medium 
the force is non zero. The total cell interface in Vl, which moves as the whole 
under the force effect, is J pdV. Hence the force acting on cell interface unit is 
EV ' pj p. Supposing the viscous law for the cell interface movement, we find a 
cell interface flux j = pv oc DV p. We see that the cells make a collective motion 
in the direction V p. This is an uphill diffusion" of the cell medium. Taking into 
account the annihilation law dp/t oc —Dp 3 , which follows from p(t) oc i -1 / 2 , we 
get g§ 

dp/dt = -71 .Dp 3 - 72 -DA/?. (5) 

Here 71 and 72 are dimensionless constants, depending on the scaling distri- 
bution function ip(z), only. Equation (5) has the form of a diffusion-reaction 
equation with the cell interface density p playing the role of concentration and 
with a negative "diffusion" coefficient. We see that, the amplitude of any small 
spatial fluctuation Sp(x, t) with wavelength A < £, where £ = lir^J 72/37ip~ 1 is 
the correlation length of the fluctuations, grows without bound. It follows that 
competing stochastic processes of cell disappearance and cell medium "uphill 
diffusion" can give rise to the formation of large-scale domains of the size A < £ 
having the density p(x, t) greater or smaller than the mean value p(t) (Fig. 2 
p3[). It also means that space derivatives of higher order of p must be taken 
into account in right side of (5) for detailed analysis of the system stability and 
more accurate estimation of £. 

Whereas the physical ideas explaining the possible macroscopic self-organi- 
sation effects are simple, the consecutive derivation of the mesoscopic evolution 
equation (5) and the calculation of kinetic coefficients 71 and 72, defining the 



effect quantitatively, turns out to be a difficult statistical task 14 1. Therefore 
the "technique" of the problem solution asks for some commentaries. Let us 
consider an ideal g clS clS el well-known example for comparison. The consistent 
description of its evolution include three classic scale levels of consideration: 
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Figure 2: Late-time configurations in Potts model of grain growth simulation. 

micro-, meso- and macroscopic ]5|, |6|, [l^, On the microscopic level the evo- 
lution is effectively described by methods of physical kinetics, i.e. by Boltzmann 
kinetic equation || ^] for the system probability distribution function / in a 
phase space, taking into account interaction and movement laws of atoms. Hav- 
ing the equilibrium solution fo for the homogeneous gas (Maxwell-Boltzmann 
(MB) distribution), it is possible to move to the mesoscopic level of the system 
description by introducing mesoscopic quantities (density, pressure, etc.) as mo- 
ments of the distribution function. In the case of small spatial gradients of the 
mesoscopic quantities, integration of kinetic equation by Enskog method gives a 
hydrodynamic equation || [l5| . If the hydrodynamic equation has steady-state 
nonhomogeneous solutions, it is possible to talk about a macroscopic description 
of the system by means of spatial patterns of size £ much greater than the meso- 
scopic size. To study the patterns formation it is convenient to transform the 
hydrodynamic equation into Langevin stochastic equation |T(| for an order pa- 
rameter. As we see above LSH distribution (4) plays the same outstanding role 
for the coarsening cellular structures in steady-state, as MB distribution plays 
for ideal gases in equilibrium state H g, 0, §. All difference between MB and 
LSH distributions is that the first is determined by thermofluctuations related 
to collisions of gas atoms and the second is determined by stochastofluctuations 
related to the process of cells annihilation. Hence the main methodical idea of 
this paper is to use LSH distribution as a basis for the multilevel description of 
cellular structure, like MB distribution is used in order to go from the kinetic 
equation to the hydrodynamic one and then to the stochastic one. This idea 
has interesting applications for other physical systems, where the scaling regime 
is reached. For example, its application allows to generalise Lifshitz-Slyozov 
theory of homogeneous coalescence for the nonhomogeneous case jl7j , to de- 
scribe the relaxation of nonhomogeneous dislocation structure during annealing 
|§ etc. 
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2 Kinetic theory (microscopic level) 



Let us consider as in Sec. 1 a nonhomogeneous cellular structure with the in- 
terface density p(x, t), weakly varying in space. We divide the volume V of the 
system into small coarse-graining elements Vl with L « A, but still includ- 
ing many cells. In the first approximation the cell dynamics of the structure 
in each Vl will be defined by the dynamics equation (1), with the curvature k 
becoming a local caracteristic - an interface curvature self- consistent mean field 
k(x, t) . The field k(x, t) in Vl at point x equals to the self-consistent mean 
of cell curvature n(t) determined by (l)-(4) for the cell structure in Vl if we 
approximately consider it as homogeneous. The field /c(x, t) is a weakly varying 
function in space as p(x, t) is. In the second approximation the dynamics equa- 
tion (1) will contain gradients of the field k(x, t). It is easy to establish this fact 
if we analyze an one-dimensional cellular structure (d = I). Let us cell interfaces 
are located at the points Xi(i — 0, ±1,±2, . . .) of x-axe. We have the cell radii 
a.W = (xi+i — Xi)/2 , the cell curvatures — 1/a^ and the cell centres coordi- 
nates a;W = (x i+ i+Xi)/2. Each interface is acted upon by attraction of adjacent 
ones (the interface tension of the cells) according to the law Fi = J?(rW — k^ -1 )) 
with the result that the cell changes its position dx^ jdt = D(k^ +1 ^ — k^ 1 ')/2 
and its size da^ jdt = D(k/- 1+1 ^ + k^~ x > — 2«W)/2. If we now introduce a cur- 
vature exact field K e {x, t) as a step function, which is «W if x, L < x < CEi+i , we 



can rewrite dx^/dt = D[n e (x^ - a&,t) - n e {x^ + a^,t)]/2 and da^/dt = 
D[n e (x (l ^ ,t)+K e (x^ - , t) - 2/a (i) ]/ 2 - Replacing the exact field K e (x, t) 
by the mean field n{x,t) and taking into account that a << L << A we get 
respectively the cell velocity and the cell size change rate 



It is easy to understand that the abovementioned definition of the curvature 
mean field fc e (x, t) holds for nonhomogeneous cellular structures in two (d = 2) 
and three dimensions ( d = 3). Expressions (6) and (7) are also valid. Really, 
each cell of the structure may be approximately thought of as having the form of 
sphere. Then the all cellular structure configuration may be described through 
the cell centre coordinates and the cell radii. Let us choose a cell of radius a with 
centre in a point x and an interface point with coordinate x' = x + a(x', t) of the 
cell. The velocity w(x', t) of the interface point affected by the interface tension 
force F(x',i) = — E[l/a(t) — K e (x', i)]a(x', t)/a(t), is determined by the viscous 
law w(x',i) = MF(x',i)]. Replacing K e (x, t) by k(x, t) we can approximately 
calculate the cell velocity dx/dt = § w(x.' ,t)ds' / (ildda d ^ 1 ) « aWw/d and the 
cell size change rate da/ 'dt = § u>(x', t)ds' /(fldda d ^ 1 ) ps w + a 2 Aw/ (2d) and 
find expresions (6) and (7). 

We see that the cells move in a coordinate- size space {x, a} in the process 
of coarsening. These movements are conveniently described by kinetic equation 



for the structure probability distribution function /(x, a,t) in the cell coordi- 
nate-size space with normalisation condition 



dx/dt = DaX7n/d, 
da/dt = -D(l/a - k) + Da 2 AK/(2d). 



(0) 
(7) 



df/dt + d(fdx/dt)/dx + d{f da/dt) /da = 



(8) 




(9) 
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Equations (8) and (9) have the physical sense of continuity equation in the 
coordinate-size space and of conservation law of the total volume of all cells of 
the system. They are natural generalisations of (2) and (3). For the macroscop- 
ically homogeneous cell structure the set of equations (6)- (9) comes to (l)-(3) 
and has the solution (4) with <p(z) = [dz/(2 - z) d+2 ](2e) d exp[2d/ (z - 2)] for 
< z < 2 and ip(z) = for z > 2 ; /i, = J z l ip(z)dz . The mean structure 
variables obey simple scaling laws: a(t) = fii^/Dt; n(t) = ((14-2/ (J>d-i) / vDi", 
n(t) = l/(n df i d (VDi) d ) and p(t) = d^ 1 /(2 f x d VDt). 

To find /i(x, a, t) — fo(a, t) + Sf(x, a, t) in the first-order approximation for 
small parameters | Vp | / p 2 « 1, it is necessary, as in well-known Enskog 
method [0], to substitute /i(x, a,t) in (6)-(9) and linearize them, taking into 
account that <5/(x, a,t) << f (a 1 t). I hope in a future work to realise this 



possibility, but this is proving at present to be a difficult task 14 



3 Hydrodynamic theory (mesoscopic level) 

Equation (8) provides the kinetic (microscopic) description of the cellular struc- 
ture coarsening. To get the hydrodynamic equation let us rewrite (4) as 

/o(o,t) =ap d+1 <p(f3ap), (10) 

where a — (p d ) d /[fl d (dp d -i) d+1 ] and {3 = p d /(dp d -i) are constants, and use 
again as in Sec. 1 and 2 the standard coarse-grain averaging procedure Q. 
Namely, let us introduce a mesoscopic coarse-graining length L{t) » d(t), 
so that a mesoscopic coarse-graining volume Vl = Q d L d includes many cells. 
Let us divide the volume V of cellular structure into elements of volume Vl- 
For small gradients of mesoscopic variables k(x, t), <z(x, t), p(x, t), and n(x, t) we 
may suppose that a local homogeneous steady-state is reached in each separate 
element Vl after a time tl = L 2 / D much greater than the cell annihilation time 
T a = a 2 /D. For all this, the cellular structure as the whole (macroscopically) 
is not in the homogeneous state. Then the distribution function /o(x, a, t) in 
Vl at point x can be assumed to be a local steady-state function, equal to LSH 
distribution (10) for the homogeneous cellular structure with density p(x, t), 
that prevails at the point x . This supposition is completely analogous to the 
local system equilibrium one in linear thermodynamics ^) . Multiplying (8) by 
Sldda 11-1 /2, (6) by Vl d da d ~ 1 fo/2, (8) by VL d a d , integrating them by a and using 
conservation law (9), we make the averaging procedure. After sufficiently heavy 
calculations we obtain, respectively, the equation of cell medium continuity , an 
analogue of Darcy 's law in hydrodynamics and an equation similar to the state 
equation in gas dynamics |15| : 

dp/dt + V(pv) = -D[ ai p 3 - a 2 p 2 n - a 3 A/t], (11) 

v = (D/2)p~ 1 X7n, (12) 

P 3 k = l3 lP 4 - I3 2P Ak + foVnVp. (13) 

The set of three hydrodynamic equations (11)-(13), containing three hydrody- 
namic quantities of the cell medium: the density p(x, t) , the velocity v(x, t) — 
(Sldd/p) J Q a d_1 (dx/di)/da, and the curvature K(x,t), forms a closed set. The 
role of the quantities p(x, t) , k(x, t) and v(x, t) are completely analogous to the 
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role of gas density, pressure and velocity in viscous gas dynamics. If initial and 
boundary conditions are known, the cell medium movement can be determined 
by traditional hydro-gas dynamics methods p"5| for any continuous distribution 
of p(x, i) , n(x,t) and v(x, t). By substituting (12)-(13) in (11), the set can be 
reduced to evolution equation (5). In (5), (11)-(13) ai, «2, fli, /?2> 7i, 72 are cal- 
culated constants. For d = 2: ct\ — 1.54, a® — 1.11, a 3 — 0.25, j3\ = 1.11, /3a 
0. /3 3 = 0.267, 7i = 0.308, 72 = 0.277. For d = 3: ai = 1.16, a 2 = 1.44, a 3 = 0.5, 
/3i = 0.718, & = 0.132, /3 3 = 0.263, 71 = 0.129, 72 = 0.136. Full expressions are 
sufficiently bulky (e.g. ct\ = 4(d — l)p 2 d pd--i,l {d 2 p\_i)) and will be published 
elsewhere fT^f . 

4 St ochasto dynamic theory (macroscopic level) 

To analyse the dynamics of macroscopic spatial perturbations of the cell me- 
dium, let us introduce a relative value of deviation of the cell interface density 
from its mean macroscopic value ip(x,t) — 8p(x,t)/p(t) << 1 , where p{t) is 
the uniform solution of (5). Substituting p(x.,t) = p(t) + 6p(x,t) in (5) and 
coming to new time and space variables dt = 3-fiDp 2 (t)dt = — d In p 3 (t) and 
<9x = p(t)dx (intrinsic system variables 0), we obtain 

dip/dt=-5F/5ip (14) 

with T{ip} = J dk[Aip 2 +C(VV>) 2 ]/2 , A = 1 and C = -72/(371) (C = -0.300 
for d = 2 and C = —0.351 for d = 3 ). Taking ip(5c,i) in the form ip(x,i) = 
?A(k, a>)e*( kx-t ' i '), we get the dispersion relation uj(k) = — icr(k) , where 

<r(k)=.4 + Cfc 2 (15) 

is the damping coefficient. Including in ^{ip} the term — /i?/> describing the effect 
of an external field h(5t, t) find the generalised susceptibility 

x(k,w) = i/[?(k)-<w]. 

As all mesoscopic quantities, ^ is a fluctuating quantity. But on writing 
(14) we have neglected these fluctuations. The quantity ip is a sum of a number 
of microscopic random quantities. Each of them corresponds to the contri- 
bution of a cell being in Vl- In view of central limit theorem of statistics, 
ip has Gaussian distribution function with the mean value < ip >l= and 
the dispersion < if; 2 >l= v/\VlN] , where v = P2(d-i)/ Pd-i — 1 = 0.109 
for d = 2 and v = 0.373 for d = 3) @. Equation (14), after taking into 
account the fluctuations, becomes Langevin stochastic equation with 
being a stochastodynamic potential [ p^| . The noise £ is Gaussian noise, with 
< C >= and the correlation function < £(x, i)((5t', t') >= 29S(x - x.')S(t - ?) 
, defined by the fluctuations at scales a « L and r Q << tl. Here = 
fi d (d/2) d - 1 (^_ 1 /^- 1 )( M2(d _i ) / At 2_ 1 -1) (fl = 0.619 for d = 2 and = 8.09 for 
d = 3 ) The forms of evolution equation (14) and of the noise £ guarantee 
that the steady-state probability distribution functional is given by Boltzmann 
formula P{i>} oc exp[— T {ip} / 9] |^9|. We see that our cell system is rated in 
the class of so-called potential systems with potential ^F{4>} having Landau- 
Ginzburg form |^|, ^ |l6| and the quantity ip plays a role of an order parameter. 
It is not surprizing because in general sense the system is in perfect analogy to 
a thermodynamic system with a nonconserved order parameter ip. So the forms 
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of evolution equation (14) and of potential T{^p} are dictated by simple symme- 
try macroscopic considerations. Instead of the thermal noise in thermodynamic 
systems we have here the stochastic noise, created by the continuous process of 
cells annihilation. In virtue of reasoning, fluctuation-dissipation theorem may 
be applied, according to which the generalised susceptibility %(k, Q) determines 
the spectral correlation function of fluctuations as S(k,£j) = (20/<D)ImxQc, uo) . 
The structure function has usual Ornstein-Zernike form S(k) = 9/a(k) H ||. 

From (15) we see, that the order parameter modes i/}0s., t) with k > d—AjC 
are unstable. This is every indication that highest order terms of the expan- 
sion of 3-{4>} in powers of V?/> must be considered for detailed analysis of the 
system stability. The derivation of the exact expression for these terms by 
Enskog method is a very heavy task || [l4| not realized to the present. We 
proceed now differently. Let us use the fact that our system is potential with 
the nonconserved order parameter ip . Then we may consider as some 

phenomenological potential in the abovemensioned Landau sense. This poten- 
tial have necessarily the general form — J dx[Aip 2 + C(Vijj) 2 + T>(Aip) 2 ]/2 
with P > |, |1| ^0|. The phenomenological coefficients A,C and T> can be 
found from experiments, computer simulations or by a direct calculation as we 
have made above for A and C . As a consequence we get the damping coefficient 
in the form 

or(k) = A + Ck 2 + Vk 4 . (16) 

We do not know the real value of V , but we are able to forecast the possible 
scenarios of system evolution Q. If A > 0, C < and V > C 2 /(4A) the 
evolution equation (14) has the single stable solution i/jq(x) = , and modes 
ip(k, t) with any k decay with time. Therefore the system is stable with respect 
to large-scale fluctuations of ip(x,t) near i[>o(x.) = . The modes ip(i.,t) with 
k = kc = y/—C/(2T>) = 27r/£ , where £ is an intrinsic correlation length, have 
the lowest decay rate. That is why under the action of the stochastic internal 
noise C > created by a continuous process of cells annihilation, space-correlated 
fluctuations arise in the system. The system is in a modulated state (short range 
ordering) and the structure function 

5(k) = 0/{A + C~k 2 + v¥). (17) 

has a maximum at k = kc (Fig. 3), reflecting the abovestated fact that the 
modes ip(ke,i) decay with the lowest rate. If A > 0, C < and < V < C 2 /{AA) 
, the evolution equation (14) has the unstable solution ^o(x) = and stationary 
one-, two-, ... d -dimensional space-periodical solutions with wavelength A = £ 
[p0[ . Under the action of the stochastic internal noise ( the system comes from 
the homogeneous state into one of these patterned state (long range ordering). 
The question of existence of the first or of the second type of ordering depends 
of the value of coefficient T> . 



5 Discussion 

We see that in the process of structure coarsening the macroscopic cell interface 
density p(t) decreases with time as p(t) oc i -1 / 2 and the mesoscopic density 
p(x, t) = p(t)[l + 6p(x,t)/p(t)] also decreases (Fig. 1). The competition of two 
stochastic processes of the cells annihilation and of the cells collective motion 
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Figure 3: Scaled structure function of fluctuations S(k)/S(0) = l/[l + (C /A)k 2 + 
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creates space-correlated fluctuations of the cell interface density. These fluctua- 
tions are the adjacent domains of size £ with p(x, t) greater or smaller than p(t) 
(Fig. 1,2). In other words, they are the domains in which the cells have the 
mean size a(x, t) greater or smaller than the macroscopic mean value a(t) . The 
amplitude of relative fluctuations <5p(x, t) / p(t) and the domain size (characteris- 
tic wavelength of fluctuations), divided by the mean cell size £,/a(t) oc £p(t) = £ 
stay constant. These domains can be short-range ordered (Fig. 1) or long-range 
(Fig. 2) ordered in space. The short-range ordering of cell interface density 
fluctuations resemble the fluid density fluctuations ordering in the well-known 
phenomenon of critical opalescence. The long-range ordering of cell interface 
density fluctuations resemble to the effect of spinodal decomposition in thermo- 
dynamic systems. 

This is just a moment to return to the question discussed in Sec. 1 and 
2 about basing of the self-consistent mean field approximation for description 
of the cellular systems. The condition of its validity is the condition that the 
relative value of field fluctuations <5k(x, t)/n{t) m Sp(x, t)/p(t) « V( x i t) is small. 
More priscisely, the fluctuations averaged over the coarse-graining volume Vl 
must be small:< ip 2 (x,t) >l~ v/[VlN(tc, t)] « 1 . Choosing the volume 
of a first coordination sphere as a minimum possible value of Vl and taking 
into consideration that n — VlN is the number of cells in Vl , we see that 
this condition is well satisfied both for d = 2(y = 0.109, n w 7) and d = 
i(i> = 0.373, n fa 13). The general condition of validity of mean-field theory 
for the systems with Landau-Ginzburg potential is well-known Ginzburg 

criterion: fluctuations of ip averaged over the correlation volume must be small 
< ip 2 (x.,t) >j= v/[V£N(x,t)] « 1 . This criterion is obviously fulfilled if 
L < £ . Uniting this condition with the foregoing condition of existence of 
stochastodynamic potential , we have finally a « L < £ . This is just 

the familiar condition of the choice of the coarse graining size L , which is also 
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the condition of validity of the self-consistent mean field approximation in phase 
transformation theory. It can be satisfied if a < < L < £ . In coarsening cellular 
systems this is always the case for the short range ordering (D > C 2 /(4A)) 
and this is so indeed if T> >> —C/(8ir 2 ) for the long-range ordering (0 < T> < 
C 2 /(AA) ). 

In spite of that the ordering effects of cell interface density fluctuations 
can be evidently established by visual analysis of different cellular structures 
in experiments (Fig. 1) and computer simulations (Fig. 2), I am not aware 
of publications, where the structure function (17), describing the effect quan- 
titatively, has been directly obtained. I feel that it is related to two reasons. 
Firstly, it must be noted that hydrodynamic and stochastodynamic theories de- 
veloped here are valid in the long wavelength limit k = kp{t) < 1 , when we can 
guarantee the possibility of application of the cell medium conception and the 
self-consistent mean field approximation. In the long wavelength limit k < 1 the 
maximum of SQs) is hardly observed: if fee << 1 the maximum is small (Fig. 
3); if kc ~ 1 it is hardly separated from the other maximum ko rs 2 correspond- 
ing to the trivial effect that the mean distance between neighbouring interfaces 
is 2a . This later maximum cannot be described within our long wavelength 
approach and is not shown in Fig. 3. Therefore the more strong confirmation 
of existence of the large-scale fluctuations of cell interface density is the obser- 
vation of the negative curvature of S(k) in the macroscopic k — > limit (Fig. 
3). However it demands of a laborious statistical analysis or a scattering study 
of very large cellular systems. In actual practice the finite-size induced effects 



stays too essential 21 . Secondly, in computer experiments, using usually " 
q- states Potts model" [|l0|, Ol, the structure function S'(k) is calculated as a 
Fourier transform of the correlation function of a phases volume density. If all 
phases have the same equilibrium energy, and the boundaries separating the 
phase domains, are sufficiently distinct, the system state is completely defined 
by the total interfacial area. In consequence, the cell interface density p(x, t) 
and the corresponding structure function (17) is a much more sensible charac- 
teristic of the system than the phases volume density and the corresponding to 
it structure function. Quite recently [ ^2| it has been demonstrated by computer 
simulation of normal grain growth that a size correlation between neighbouring 
grains exist in grain structure. Unfortunately the spectral correlation analy- 
sis of the cell interface density fluctuations by means of the structure function 
(17) has not been made and direct quantitative comparison of these results with 
present theory results is difficult. Still there is no doubt that we are dealing with 
the first observation of the short-range ordering of cell interface density fluctu- 
ations in cellular structures. New experiments or computer simulations with 
careful space spectral analysis of cell interface density fluctuations are needed 
to confirm the existence of dynamic patterning in cellular structures. 



6 Summary 

In summary, the problems of multilevel description of cellular structure evolu- 
tion are considered from positions of kinetic, hydrodynamic and stochastody- 
namic theories of nonequilibrium statistical systems. It is demonstrated that 
the effective method to solve the problem is to describe the cellular structure on 
the microscopic level by the probability distribution function in the coordinatc- 
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size space, i.e. in the space of the cell coordinates and the cell sizes. The 
kinetic equation for the weakly nonhomogencous cellular structure is written 
and the development to hydrodynamic equations for a continuous elastic cellu- 
lar medium with viscous motion law is realised. These non-linear differential 
equations contain three mesoscopic quantities of the cell medium (the cell inter- 
face density, the cell medium velocity and the cell interface curvature) and are 
analogous to the continuity equation, Darcy's law and the equation of state in 
gas dynamics. For continuous flow of the cell medium the hydrodynamic equa- 
tion set is reduced to a single evolution equation for the cell interface density. 
It has the form of a diffusion- reaction equation with a negative " diffusion" co- 
efficient and with the cell interface density playing the role of " concentration" . 
It is shown that the relative value of deviation of the cell interface density from 
its mean macroscopic value is a suitable order parameter of the system. The 
order parameter dynamics are described by Langevin stochastic equation with a 
stochastodynamic potential having Ginzburg-Landau form and with a stochas- 
tic noise, created by the continuous process of cells annihilation. Analysis of 
the order parameter dynamics in coarsening reveals a new effect of macroscopic 
patterning in the cell medium: competition between two stochastic processes of 
the small cells " eating up" by large ones and of the " uphill cells diffusion" leads 
to the creation of space-correlated fluctuations of the cell interface density. 
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